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ABSTRACT 

The following research results are based on develop- 
ment of an approach previously proposed by the authors 
for optimum nozzle design to obtain maximum thrust. 
The design was denoted a Telescope nozzle. A Telescope 
nozzle contains one or several internal designs of certain 
location, which are inserted at certain locations into a 
divergent conical or planar main nozzle near its exit. 
Such a design provides additional thrust augmentation 
over 20% by comparison with the optimum single nozzle 
of equivalent lateral area. What is more, recent experi- 
mental acoustic tests have discovered an essential noise 
reduction due to Telescope nozzles application. In this 
paper, some additional theoretical results are presented 
for Telescope nozzles and a similar approach is applied 
for aeroperformance improvment of a supersonic inlet. 
In addition, a classic gas dynamics problem of a simi- 
lar supersonic flow into a plate has been analyzed. In 
some particular cases, new exact analytical solutions 
are obtained for a flow into a wedge with an oblique 
shock wave. Numerical simulations were conducted for 
supersonuc flow into a divergent portion of a 2D or ax- 
isymmetric nozzle with several plane or conuical designs 
as well as into a 2D or axisymmetric supersonic inlet 
with a forebody. The 1st order Kryko-Godunov march- 
ing numerical scheme for inviscid supersonic flows was 
used. Several cases were tested using the NASA CFL3d 
code based on full Navier-Stokes equations. Numerical 
simulation results have confirmed essential benefits of 
Telescope design applications in propulsion systems. 

I. INTRODUCTION 

Several well-known experimental results show essen- 
tial acoustic benefits in the application of some untra- 
ditional nozzle designs. For example, nozzles with rect- 
angular or elliptic cross section in the supersonic part 
produce less jet noise than round nozzles designed for a 
fixed Mach number at the nozzle exit ( i.e. with uniform 
flow at the exit and pressure coinciding with the flight 
static pressure outside the exhausting jet). Thus, the 
theoretical perfectly shock free jets are ’’noisier” than 
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at least partially underexpanded (or overexpanded) jets 
with possible internal shocks. Moreover, the experimen- 
tal research of Ahuja, Krothapalli et al. [1,2] has shown 
that inserting disturbing elements into supersonic jet 
flow: slots, finger, tabs etc., can reduce jet noise (and 
screech tones) in spite of the presence numerous strong 
and weak shock waves. This contradicts the traditional 
view on the considered phenomenon. A reasonable ex- 
planation for these facts would be the appearance of 
more effective mixing and destruction of the regular 
cell-shock structure in the weakly underexpanded jet. 
Inside such a jet, the weak barrel-shaped shock waves 
are always present and these shock waves are the main 
sources of the oscillatory processes in the jet. In the reg- 
ular almost parallel co- annular mixing layers, unstable 
longitudinal waves are excited, and noise is produced in 
the fixed direction of the jet axis ~ 145°. Of course, the 
presence of shock waves in the jet exhaust, especially for 
a supersonic nozzle, can lead to some dangerous side ef- 
fects and performance penalties. 

Developing previous ideas for jet noise reduction, two 
novel concepts were proposed in the papers [3-6]. This 
first concept is denoted as the Telescope nozzle, for 
it consists of several internal nozzle surfaces that are 
arranged in a telescope fashion. The second concept is 
denoted as the Bluebell nozzle, based on the flower- 
like shape of its external jet plume. Bluebell nozzles 
utilize both chevrons and corrugation in its nozzle ge- 
ometry. Each concept is capable of achieving a thrust 
performance greater than the standard baseline conic 
or 2D plane convergent and convergent- divergent (CD) 
nozzles. The improved performance for Bluebell noz- 
zles occurs due to increase in nozzle internal surface 
area while maintaining nozzle-projected area equiva- 
lent to the baseline reference nozzle. Small scale and 
large scale acoustic tests of different modifications of 
Bluebell nozzles were conducted at the NASA Langley 
Research Center and Central AeroHydrodynamics In- 
stitute (TsAGI) in Moscow, Russia. These tests have 
shown essential acoustic benefits Bluebell design ap- 
plication in supersonic regimes as well as in subsonic 
regims. For example, the experimental tests of several 
Bluebell nozzle designs ([3]) have shown noise reduction 
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relative to a CD round nozzle with design exhaust Mach 
number M e = 1.5. The best design provides an acoustic 
benefit near 4dB with about 1'/, thrust augmentation. 
Below, we consider only the first (Telescope) concept 
with the goal of the Telescope nozzle design optimiza- 
tion for the muximum nozzle thrust, with the intent 
application of this concept to propulsion systems, es- 
pecially, for a supersonic engine inlet. Detailed infor- 
mation about the second (Bluebell) concept is in the 
papers [3-6] and in the patent [7]. 

II. PLANE ELEMENT IN SUPERSONIC FLOW 


2.1. The thrust on a plate element with an 
oblique shock wave and Prandtl-Meyer rarefac- 
tion flow. A divergent flow can act on a plate or airfoil 
inserted into a flow so that a resulting force is directed 
against the flow. This effect is used for thrust by su- 
personic nozzles. Conversely, a uniform flow produces 
only drag for bodies and airfoils. Inserting a conical or 
wedge-shaped nozzle inside the divergent part of an ex- 
ternal nozzle so that the integral of the pressure on the 
low side of the inserted surface is greater than on the up- 
per side produces increased thrust. There is an optimal 
angle of the plate that provides the maximum thrust at 
each point of a divergent flow. The most efficient inter- 
nal design is produced from a pattern that looks like a 
telescope with extending tubes. The optimal number of 
internal designs is defined through dependence on the 
Mach number at the nozzle exit, M e . Telescoping de- 
signs must be located so that the compressible waves 
formed by interaction of a flow with this design would 
be passed on to the upper side of the next lower tele- 
scoping part. The best result will be produced by such 
a set if the external design inclination increases down- 
stream. Computations show that a significant thrust 
benefit from the Telescope nozzle occurs with an exter- 
nal telescoping design, using either wedge, conical or 
optimal contour shapes, and also in the case of a plug 
application. 


This effect can be demonstrated by consideration of 
the classical steady supersonic flow which forms when 
an infinite uniform supersonic flow diverges at the de- 
viated wall (convex angle wall). The schematic wave 
picture of such a flow is shown in Figure 1. A super- 
sonic flow with Mach number, M = >1 is flowing 

along a rectilinear wall EO. After point O, the wall is 
turned at some angle 6 0 . The flow velocity increases in 
the centre wave so that all flow parameters depend only 
on the angle p. All flow parameters can be calculated 
by: 

A 2 = 1H Ar sin 2 (atp), — = (1 - a 2 A 2 ) 6 (1) 
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w r = -sin(a<p), w v = cos(aip) (3) 

sina = 6 = a + <p - ^ (4) 

where 

2 1 
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Here the angle p is measured from the straight line <p = 
0 which is inclined to the initial wall direction by the 
angle 9?^,; <5 is the angle between local velocity and initial 
wall direction. 

Let a plate element be inserted into this Prandtl- 
Meyer rarefaction wave. Then depending on this ele- 
ment’s location and inclination to the local velocity vec- 
tor, W (wriWy), two different cases are possible. The 
first, when the element is located in the uniform flow 
region (the region 5 in Figure 1). Then the flow at the 
plate is similar and the thrust (or drag) produced by 
this element can be calculated by relationships in the 
Prandtl-Meyer rarefaction wave (1-6) for one side of the 
plate and by relationships on the shock wave (BC) for 
the other side: 


/(/?)/ tan 

WO?) ’ 


/(/?) = ^-(sinV-M- 2 ) (7) 


where /? and 7 are local angles between the flow di- 
rection (the vector W) and the shock wave and plate 
respectively. This case is analyzed in paper [3]. It 
was shown that the optimal inclination of the plate de- 
pended on the local Mach number and flow direction. 
But these results were not connected with any fixed 
flow, so that the two given parameters, Mach number 
and angle of flow direction are independent. In the case 
of rarefaction flow, these parameters are dependent by 
relationships (1-6). Calculations show that maximal 
thrust is produced by the unit plate element when this 
plate is inclined to the initial wall direction of the angle, 
7 ~ 10 °. 

In the second case the plate is located in the rar- 
efaction wave and the flow at the plate is not similar. 
The unknown oblique shock wave is not rectilinear, and 
the flow behind it is not uniform and parallel to the 
plate direction. In some approximations, this problem 
has an analytical solution. For example, the flow on 
both sides of the plate can be calculated in acoustic ap- 
proximation for weak waves using an analytical or any 
numerical method for 2D flow. 

2.2. Oblique shock relationship analysis. Some 
theoretical analysis of shock wave relationship (7) was 
conducted with the goal to simplify and increase the 
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efficiency of the numerical algorithms. Let a plate (or 
wedge) be located at the angle, 7 to the supersonic flow 
direction. It is well known that the formula (7) defines 
this angle implicitly through use of the angle, /?, be- 
tween the shock wave and downstream flow direction. 
Usually, the plate or wedge angle, 7, is given, and the 
shock wave angle, /?, is unknown, i.e. its value is cal- 
culated by iteration of (7) or by interpolation of the 
tabled curves 7 =F(/?, Mqo) where F is some function. 
Of course, explicit formulae are preferable, especially, 
in more complicated algorithms, which calculate this 
problem solution repeatedly as, for example, in mod- 
ern numerical schemes with the Riemann solver. The 
equation (7) can be transformed to a cubic equation of 
the variable y=f (/?). For simplification, substitute also: 
x—tan 2 7, a=(/c + l)/2, m =M<^ 2 . Then this eqation be- 
comes: 

Ay 3 4- By 2 4- Cy + D = 0 (8) 

where A=a(l+x), B=m-l+(m-2a)x, C=(a-2m)x, D=mx 

Explicit solution of this equation can be given us- 
ing Cordano formulae (see [8]). In accordance with the 
analysis of the coefficients in (8), we conclude that this 
equation has three real roots. The simpliest represen- 
tation of these roots can be given by introduction of 
the auxiliary angle. Only two positive roots are in the 
range of the variable y. The greater root corresponds 
to a subsonic value of shock wave inclination, and lower 
is for a supersonic value. Since these formulae for roots 
are combersome, we omit them in this paper. 

Two important relationships can be found from (7) 
without use of its exact solution. The first is defi- 
nition of limit values for the angles 7 and 0. They 
are determined at the maximum value of the function 
0 = £(7), i.e. using the relationship: d/3/d^y = 0. 
Equating the derivative of the variable 0 on the right 
side in (7) to zero, we get the square equation for a 
variable 2 — sin 2 0: 

kz 2 — (a — 2m)z — (m 2 4 am) = 0 
and positive root is 

~ r = ^ - ™ + \J(^-m) 2 + K(m'-+am)) (9) 

so that limited shock inclination 0u m is calculated as 
fitim = asin(yfz) and the limited wedge inclination 7n* m 
is calculated using formulae (7) with 0u m on the right 
side. For hypersonic flows with Mach number, = 
oo (m=0) we have: 
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The corresponding curves for dependance of limit angles 
7 and 0 on Mach number and specific heat ratio k are 
shown in Figures 2 which were calculated using (10) and 

( 7 ). 

The second relationship is for hypersonic flows with 
Mach number, Moo = oo. In this case, the solution 
comes to a solution of the square equation relative to 
variable, z ~ sin 2 0. Designate t=tan 7 , then the roots 
of this equation are: 


* 1,2 


1 + (« + l)t ± s/i - (k 2 - 1)<3 

2(1 + f 2 ) 
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where ” plus” sign corresponds to subsonic conditions and 
” minus” to supersonic condition. The set of well known 
” apple-like” curves for dependence angles 0 = F(7) for 
different specific heat ratios are shown in Figure 3. The 
black poins with numbers 1-11 correspond to limited 
values of these angles. 


2.3. Optimum plate location. In common case, 
a plate (or airfoil) inserted into an inviscid supersonic 
flow produces a resulting force normal to the plate and 
its value and direction depend on the pressure difference 
on both sides of the plate. The nondimensional aero- 
dynamic characteristics of the plate, the thrust, T ni or 
drag, Cd , produced by this flow about the plate can be 
calculated with these four parameters: specific heat ra- 
tio K t flow Mach number, Moo, an angle a between the 
flow and thrust direction, and the angle 7 between the 
flow and the plate. Schematic geometry and designa- 
tions for such a flow are shown in Figure 4a. An angle, 
7, is measured from the upstream flow direction. If 7 >0 
and less than the limited angle 7/» m , (0< 7 < jn m ), 
the thrust (drag) is determined by the simple analyti- 
cal formulae using relationships for oblique shock waves 
and the Prandtl-Meyer rarefaction wave discussed in 
the previous sections. In this interval of the angle 7, 
for all another parameters there is an optimal value of 
the angle 7 opt , which gives the thrust maximum value. 
Aerodynamic characteristics of the unit plane element 
in supersonic flow were calculated using the created 
code for a wide range of the parameters: k, Moo, and 
a. Two examples are illustrated in Figures 4b, c by the 
curve families for the nondimensional variables (T n) y°) 
for Mach number Moo = 2 and M^ = 6 respectively. 
Numbers 1,2, 3, 4, and 5 at the curves correspond the 
angls, a = 0°, 22.5°, 45°, 67.5*, and 90°. Note, that for 
large attack angles to the thrust direction, a, maximum 
thrust is obtained at the limited angles, yn m . 

Similar results were observed for another case that 
correspons to a pure Prandtl-Meyer flow at the turn- 
ing point of the 2D nozzle wall. The schematic pic- 
ture of such flow and designations are shown in Figure 
5a. Tipical examples for Mach number, Moo = 1.5 and 
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3 are illustrated in Figure 5b, c. Here, the nondimen- 
sional thrust T n vs wall inclination angle, 7, for several 
slopes of the flow (or upstream wall) is given. Num- 
bers 1,2,3, 4, and 5 on the curves correspond to angle 
a = 0°, 10*, 20*, 30*, and 40*. Again, for small angles, 
a, there are maximum thrust values inside the interval 
G< 7 < 7 Hmt and for greater a, maximum thrust occurs 
at the limited value ju m . Similar results were observed 
with other Mach numbers. 

III. TELESCOPE NOZZLE NUMERICAL 
SIMULATION RESULTS 

The main numerical simulation results for Telescope 
nozzle based on the 1st order Krayko-Godunov march- 
ing scheme [9] for inviscid supersonic flows were des- 
cribed in the papers [3,5]. Some of them will be illus- 
trated below. In particular, it was shown that by the 
standard deformation of single axisymmetric or 2D noz- 
zle it is possible to obtain only an insignificant thrust 
augmentation, even using optimal nozzle shape, i.e. ~ 
1-4'/. by comparison with usual conical (wedge-shaped) 
nozzle. A Telescope nozzle can increase thrust signifi- 
cantly more. 

3.1 Telescope nozzle geometry. Example of two 
possible Telescope nozzle embodiments are shown in 
Figure 7a and 7b. In Figure 7a, the external nozzle 
is constructed by giving the fixed contour z=z(x) in the 
zx-plane and a cross section contour is described by the 
super-elliptical equation: 

(: yM*))" (x) + (z/*«) n(s) = 1, * = /(*) 

n(x) = 2 + H(x - x*) • n e (x - x*)/((x« - x *) (12) 

c(z) — u(r)//?(x) ~ 1 i7(x x*)-c e (x x*)/((x tf x*) 

(13) 

where the Heaviside function H(x-x,) is defined: H(x- 
x*)=0 if x 0 < x < x*, and H(x-x*)=l if x* < x < x t . 
The subscripted indicies 0,* and e correspond respec- 
tively to the nozzle inlet, throat and exit. The subsonic 
portion of the nozzle (from the inlet to the throat) has 
axisymmetric shape (a=b=l, n=2). In the supersonic 
portion (from the throat to the nozzle exit), a power n 
in (12) changes from the minimal throat value of 2 to 
the maximal exit value n e and an eccentricity c=a/b 
changes from the minimal throat value of 1 to the max- 
imal exit value c e . The nozzle contour z=f(x), in the 
plane of symmetry, y=0, is a cubic parabola in the sub- 
sonic part and then becomes rectilinear with the angle 
or = 10°. For the Telescope nozzle in Figure 6a, the 
power n(x) increases from 2 to 10 downstream from 
the throat to the exit and two plane internal designs 
located symmetrically supported by the holders into 
the external design. In Figure 6b another style of the 
Telescope nozzle is shown. The external nozzle is de- 
scribed by a super-elliptical equation with constant val- 
ues n=2, a=b=i, and the internal design is described 


by the same equation with a=l,b=0.5, i.e. cross section 
has an elliptical shape; the internal design is again sup- 
ported by holders to the external design. In the both 
cases, the nozzle contour in XZ-plane, z=f(x), contains 
rectilinear and round intervals with continuous inclina- 
tions in the points of discontinuous curvature. 

Another embodiment is illustrated in Figure 6c. The 
external design is a Chisel nozzle. This nozzle can be 
constructed on the base of any plain nozzle. For the 
simplest design in Figure 6c, dependence of radius on 
the azimuthal angle in the cross section is described by 
a periodic function r = r(<p) with a period T=2 ir/n c : 
into the first period r = r + =const for G< <p < (p\ and 
^2 < <P < T y and r = r_=const for <pi < <p < (p 2i 
where <pi = 0.5(T- A <p) and <p 2 = 0.5(T + A <p). We 
call a corrugated surface part a “cavity” or a “convex- 
ity** relatively to the internal normal to the nozzle wall. 
The cavity depth (or convexity height) defined by the 
equality Ar =r+ — r_ , increases along the nozzle center- 
line from zero at the throat, x = x,, to the maximum 
value at the exit, A r =A£f , where £ = (x_x*)/(x e — x*), 
i.e. this coincides with definition of a corrugation ampli- 
tude coefficient 6. The cavity (convexity) width, A <p= 
<p 2 — <p\ } also linearly increases (decreases) downstream 
from zero (maximum) at the throat to the maximum 
(zero) at the nozzle exit, i.e. 2<p\ = T — Ay?=Tf . For 
such a configuration, two expanded flows near the noz- 
zle wall flow into two neighbouring cavities to meet each 
other at some angle a, mutually penetrate and more ef- 
fectively mix. A flow impulse on the lateral area of the 
convexities increases the resulting nozzle thrust. 

A Chisel nozzle is very convenient to use with a Tele- 
scope nozzle as shown in Figure 8 because similar con- 
vexities allow the internal design to be maintained. In 
Figure 6c the main external design (1), a Chisel nozzle, 
is based on the cone of angle a = 10° and the internal 
design (2) has a conical surface. 

The Blubell nozzle concept application for jet noise 
reduction can be also used in the Telescope nozzle con- 
cept for external design or for internal design as well 
as for both. Such an example is shown in Figure 6d. 
There a 6-petal internal conical design is installed into 
an 8-petal external design in which three plane holders 
maintain its position. 

3.2 Numerical simulation results. The numer- 
ical simulations were conducted using a modified nu- 
merical code based on the 1st order explicit numeri- 
cal marching scheme of Kryko-Godunov [9]. Solutions 
are obtained using an arbitrary curvilinear coordinate 
system, and the marching coordinate x is chosen close 
to the local streamline. A multi-zone approach and 
nonuniform grid application were used to obtain results 
of high resolution in complicated geometric domains. 
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One 3D numerical simulation result will be illustrated 
for a Bluebell shaped Telescope nozzle. In Figure 7a-c, 
this is the Telescope nozzle with conical external design 
at an angle a 2 = 20° and of internal Bluebell design 
with frequency of eight corrugations without petals, 
6 = 0 and S 0 = 0.4. The goal of this numerical sim- 
ulation is to examine the difference between the thrust 
of the baseline round conical internal design at an angle 
= 5° and the Bluebell design thrust. Both designs 
have the same corresponding cross section areas. The 
average radius of the Bluebell nozzle is calculated by 
the baseline radius. Therefore the maximal angle of 
the Bluebell internal design in the XOZ plane exceeds 
5°, and the minimal in the other plane of symmetry 
is also less than 5°. Mach contours in the sequential 
downstream cross sections from the throat to the noz- 
zle exit in Figure 12b, c show intensive development of 
the swirling flow that is swept into both the internal 
and external flows. This will favor mixing downstream 
in the jet exhaust and will reduce jet noise. The thrust 
augmentation of the Telescope nozzle application in this 
case is ~ 10’/.. 

The thrust produced by the conical nozzle divergent 
portion as a function of the nondimensional length of 
this portion L/r* is shown in Figure 8. In this numer- 
ical simulation, the maximal nozzle radius at the exit, 
r fi , was fixed at r t = 1.75 r, where the subscript * cor- 
responds to the nozzle throat cross section, and r=r(x) 
is the nozzle contour equation in the meridional plane. 

The nondimensional thrust is given T n = 7k/(/>*c* 2 * 
r 2 ), where p is gas density and c is sound speed. The 
subscripts are c for the single conical nozzle without any 
internal design; ”sh”-for the total thrust {T s h = T* - b 
T c ); and ”t”-for the Telescope conical nozzle with one 
internal design which has the length /*=0.3L. The inter- 
nal design is a plate with the entry coordinate x° =0.7L 
end coordinate x* = x tf , i.e. at the nozzle exit. The en- 
try radius of the plate is calculated from the condition 
that this entry point is located at the same flow stream- 
line starting at the throat xj 1 = 0, r* = 0.8, so that the 
plate angle is changed as shown in Figure 10 (a*). For 
comparison, the external nozzle angle, a c is also shown 
in this Figure. It is not at the optimal location and 
length for maximal thrust. These results illustrate that 
internal design is effective for short nozzles with rela- 
tively large angle a c . 

The benefit can be more significant for Telescope 
nozzles with several internal components. Figure 9 il- 
lustrates this effect that for three internal components. 
The thrust benefit here is ~25%. Depending on the op- 
timization conditions, this value can be even higher. 

Analogous results were obtained for 2D wedge-shaped 
Telescope nozzles. Comparison of the thrust for the 


single nozzle and Telescope nozzle with four internal 
components has shown that thrust augmentation can 
approach ~27'/.. For hypersonic nozzles, this augmen- 
tation is again greater. The working efficiency of a 
Telescope nozzle grows as Mach number increases. For 
example, the thrust augmentation produced by an in- 
ternal plate can reach ~75-100'/. if this plate would be 
located in a region of constant parameters of a pure 
Prandtl- Meyer rarefaction wave. An example of nu- 
merical simulation results for the 2D wedge-shaped noz- 
zle with four internal thin airfoils is shown in Figure 
10. Here Mach contours and four streamlines are pre- 
sented. These streamlines correspond to the internal 
airfoil stagnation points and also represent the zone 
boundaries. This picture illustrates the essential bene- 
fits of internal design applications inside Prandl-Meyer 
rarefaction wave and region of the constant parameters 
at the nozzle wall. In this case, the inlet Mach number 
is Moo =2, the angle /? = 30°, and the thrust augmen- 
tation for divergent nozzle portion is h=DT/T~75'/.. 
Note that the working efficiency of a Telescope nozzle 
grows as inlet Mach number increases. For example, for 
Moo =5, the value of h can mount to lOOThis value can 
be even larger depending on the angle /?. Experimental 
measurements of small and large scale Telescope nozzles 
are being addressed in future work. These experimental 
tests are very important to conduct because the recent 
experimental acoustic tests have discovered essential jet 
noise reduction produced by nozzles with the uniform 
plate set at the nozzle exit. Such noise benefit was pre- 
dicted in the papers [3,5]. 

IV. SUPERSONIC TELESCOPE INLET 

4.1 Supersonic inlet problems. The main pur- 
pose of a supersonic inlet is to slow down gas flow and 
to transfer it from supersonic speed to low subsonic 
speed before a chamber (compressor). Simultaneously, 
the total pressure should have minimal loss for effec- 
tive combustion in a chamber. The first investigations 
and analysis of this problem took place in the 60’s. For 
2D and axisymmetric inlets, the investigations showed 
that flow total pressure loss through a set of incline 
oblique shock waves with the last normal shock wave is 
essentially less than through a unique detached shock 
wave before the inlet. Several possible inlet flow regimes 
are shown schematically in Figures lla-h: a) with two 
shock waves at the inlet plus one external at the cowl; b) 
three plus one; c) three with detached shock wave at the 
cowl; d) continuous compressive waves along a curved 
inlet surface (forebody) with detached shock wave at 
the cowl; e) with partial external compression and in- 
ternal; f) with internal compression; g) with detached 
shock wave befoe the inlet. Calculations conducted by 
G.I. Petrov, K. Oswatisch have shown that the total 
pressure loss in an ideal inlet may be not more than 
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— 3-5% for a shock wave system with equal intensity. In 
real inlet flow, separation zones can be formed at the 
sharp change of centerbody inclination or at the point 
of interaction of the shock wave with the boundary layer 
at, the forebody or cowl. This makes the inlet charac- 
teristics worse. Detailed analysis of the supersonic inlet 
problem is in G.N. Abromovitch’s book [10]. Note that 
most of the optimization theories and numerical simula- 
tion methods for improving inlet efficiency do not take 
these effects into account. 

Separation and inlet drag are important obstacles 
for inlet efficient work. To reduce these effects, applica- 
tion of 3D corrugated surfaces similar to those that were 
used for improving nozzle designs may be employed. For 
example, a star shaped forebody or its smooth modifica- 
tion can reduce forebody drag. Also, it is known from 
hydrodynamic stability theory (Lin C.C., and others) 
that 2D velocity distributions in boundary or mixing 
layers are less stable than corresponding distributions 
in 3D cases because there is one additional degree of 
freedom for perturbation amelioration. Semi-empirical 
separation criteria show the same phenomenon. Several 
unusual curvilinear surfaces were proposed and tested 
experimentally many years ago by Russian scientists. 
However, such shapes have not been used in the avia- 
tion industry and require further research. Preliminary 
estimations are very promising. 

4.2 Optimum cowl shape analysis. The impor- 
tant inlet problem is to reduce inlet drag, i.e. forebody 
and cowl. Our calculations, using analytical relation- 
ships for shck waves have shown, that inlet drag may be 
reduced through judicious choice of forebody and cowl 
location and shapees. A schematic draft of one possible 
configuration is shown in Figure 14h, where instead of 
turning oblique shock waves around a cowl front edge, 
a forebody corner point was used as the turning point. 
The problem is optimization of the length and angle of 
cowl front door for a 2D inlet of an air-breathing propul- 
sion system. A systems of three and four oblique shock 
waves forming at the designed inlet for fixed Mach num- 
ber with a one and double wedge-shaped forebody were 
analyzed. The leading edge of a cowl door is the point of 
crossing of two or three oblique shock waves and the last 
passes the forebody corner point. The last shock wave 
joins this corner point and back edge of the cowl door. 
The shock wave system compresses and turns the initial 
supersonic flow and again deflects it to the same initial 
direction. The forebody drag, Cg, and total pressure 
losses , e = (p 0 -Poe)/p<>, in the cross section behind 
the forebody corner point were calculated. Here p 0 and 
p ot are the total pressures behind normal shock waves 
upstream of the inlet and behind the corner point. The 
results are presented as curves for A/^ and 0 (where /? 
is an angle between flow and plane cowl door) in the 


variables (Cg, ) and (Cg, M^) respectively. The main 
conclusion from this calculation is that for each Mach 
number there is an optimal value f3 opt >0 which pro- 
vides a maximum value for working efficiency param- 
eter of such 3 or 4 shocks inlet, f = (Cg The 

main conclusion for this problem is that the frontal cowl 
cap can create the additional engine thrust, and wedge 
shaped forebody drag can be reduced up to 15-20%. Nu- 
merical simulations by K-G marching code have shown 
approximately the same estimations for axisymmetric 
or 2D inlet with a conical or wedge-shaped forebody 
nose. The benefits increase with increase a cone (wedge) 
angle and increase in flight Mach number. 

4.3 Some numerical simulation results. The 
most previous untraditional nozzle designs discussed above 
cam be employed for a supersonic inlet improvement. In 
particular, a Telescope nozzle and all results of theoret- 
ical analysis of this concept are useful. In this case, the 
energy of the turned flow along the forebody wall can 
be used for creation of additional thrust as in the previ- 
ous case with a cowl door. As in the first problem, the 
mutual locations, sizes and angles of the internal plates 
(thin airfoils) are very important for efficiency of the 
application. Optimal values of geometric parameters 
were determined from multi-parametric numerical sim- 
ulations based on the modified marching K-G code. The 
effect of four thin airfoils installed at the minimal cross 
section (nearby of the corner point) is illustrated Figure 
12. Here Mach contours and corresponding streamlines 
are shown for the 2D Telescope inlet with a wedged fore- 
body. This design provides a forebody drag reduction 
of 25%. 

Some previous designs were tested numerically with 
the purpose of upstream viscous effects influence to 
the main conclusions. Partially, viscows effects were 
counted using boundary layer correction in thrust cal- 
culation. These estimations were conducted on mid- 
dle size vehicles for low flight altitudes. Two numerical 
codes were used: the NASA LaRC CFL3D code [11] 
and CRAFT-Tech code [12]. Both codes are based on 
an implicit upwind 2nd order numerical schemes (ENO 
versions) for solution of the full unsteady and steady 
Navier-Stokes and Euler equations. An example of a 
such result is shown in Figure 13. For these conditions 
(high Mach and Reynolds numbers), the tests have con- 
firmed the main conclusions of Telescope nozzles and 
inlets efficiency which were made on the basis of an- 
alytical solutions and numerical simulations using the 
simplified marching scheme Krayko-Godunov [9]. 

Obviously, the same approach is applicable for other 
designs, such as transition sections inside variable cross 
section supersonic tunnels, blunt bodies with several 
ring-shaped sheets, projectiles, etc. 
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V. CONCLUSION 

Theoretical analysis and numerical simulation results 
were obtained for nozzles and supersonic inlets with the 
goal of aeroperformance improvement. The designs in- 
vestigated are based on development of the approach 
proposed by the authors for optimum nozzle design for 
obtaining maximum thrust. Such a design was denoted 
a Telescope nozzle. A Telescope nozzle contains one or 
several internal designs at certain locations in the diver- 
gent conical or planar main design near its exit. Such 
design provides additional thrust augmentation over 20- 
307, by comparison with the optimum single nozzle of 
the equivalent lateral area. Recent experimental acous- 
tic tests have discovered essential noise reduction due 
to Telescope nozzles application as well. Some addi- 
tional theoretical results were presented for the Tele- 
scope nozzle and a similar approach was applied for 
aeroperformance improvment of a supersonic inlet. At 
the same time, the classic gas dynamics problem of a 
similar flow at the plate in a supersonic flow has been 
analyzed. In some particular cases, new exact analytical 
solutions were obtained for a flow at the wedge with an 
oblique shock wave. Numerical simulations were con- 
ducted for supersonuc flow into a divergent portion of 
2D, axisymmetric and 3D nozzles with several plane, 
conuical or corrugated designs as well as into a 2D or 
axisymmetric supersonic inlet with a forebody. The 1st 
order Kryko-Godunov marching numerical scheme for 
inviscid supersonic flows was used with boundary layer 
correction in thrust calculation formulae. Several cases 
were tested using the NASA CFL3d and Tech-CRAFT 
codes based on the full Navier-Stokes equation. Numer- 
ical simulation results have confirmed essential benefits 
of Telescope design applications in propulsion systems. 
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Figure 1. A similar flow in a Prandtl-Meyer rar- 
efaction wave; the thrust of the unit plate element 
(BD) in this wave. The regions and nomenclature: 1)- 
supersonic inviscid uniform flow at the rectilinear wall, 
2)- steady rarefaction wave, p-angle of the 1st family 
characteristic, a- local angle between this characteris- 
tic and a streamline; 3)-uniform flow at the plate be- 
hind' the oblique shock wave (BC), 0- angle between 
the plate element and shock wave; 4)-Prandtl-Meyer 
rarefaction flow at the plate element, 7-angle between 
this element and the streamline; 5)-uniform flow at the 
wall behind its turn, 6 0 -turn angle. 



Figure 3. Analytical solution for supersonic flow at 
the wedge with oblique shock wave for Mach number, 
Mas — oc. Shock wave angle, 0 vs wedge angle, a. 




Figure 2. Analytical solution for supersonic flow at the 
wedge with oblique shock wave: a)-Shock wave angle, 
0, vs Mach number function, for different specific 
heat ratios: /c=1.0, 1.1, ...,1.9, 2.0; b- Wedge angle, 7, vs 
A/“ s ; c)-Shock wave angle, 0, vs wedge angle, 7. 
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Figure 4. Unit plane element thrust, T n , vs flow 
angle a, plate slope, 7, and flow Mach number, A/^; 
Analytical solution for supersonic flow at the plate with 
oblique shock wave and Prandtl-Meyer rarefaction wave: 
a)- Draft and designations; b) and c)- nondimensional 
thrust vs plate slope, 7 for ^=2 and 6 respectively. 
Numbers: l, 2, 3 ,4, and 5 are for flow attack angle a=0 0 
22.5, 45.0, 67.5, and 90.0 (deg.). 


Figure 5. Unit plane element thrust, T n , vs flow 
angle 7co » plate slope, a, and flow Mach number, Moo\ 
Analytical solution for supersonic flow at the turning 
wall with the Prandtl-Meyer rarefaction wave: a)- Draft 
and designations; b) and c)- nondimensional thrust vs 
plate slope, a for A/00 = 1.5 and 3 respectively. Numbers: 
l, 2,3, 4, and 5 are for flow attack angle 7oo=0, 10, 20, 
30, ami 40 (deg.). 
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Figure 6 The Telescope nozzle embodiments: a- with Figure 7. The Bluebell-Telescope nozzle with 8 cor- 

the rectangular cross sections; b-with the elliptic cross rugations of the internal design without petals. Mach 

section, c)- Chisel-Telescope nozzle; d)-Bluebell Tele- contours: a)-inthe plane of symmetry XOZ; b)-in cross 

scope nozzle; The internal designs (number 2) are sup- section x=2.5, and c)- in cross section at the nozzle exit, 

ported into the external designs (number 1) by the hold- x=3.28. 

ers ( number 3). 
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Figure 8 Conical Telescope nozzle: thrust, area, an- 
gles vs length divergent portion: r e 1.75, AL=0.3L. Sub- 
scripts: ”c”-for single conical nozzle; ”sK”-total thrust; 
”t”-for Telescope nozzle with one internal design. Sym- 
bols: T-thrust; a-half-cone angle; S-lateral area. 


Figure 11. Supersonic inlet schemes. 
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Figure 9. Conical Telescope nozzle with three internal 
plane designs at the nozzle exit. Mach contours. 


Figure 12. 2D Telescope inlet with the four internal 
designs. 
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2D INLET NUMERICAL SIMULATION 
NAVIER- STOKES NASA CFL3D CODE 
L (2D Case, M =2, Re =10 5 ) 
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Figure 10. 2D Telescope nozzle with four internal de- 
signs at the exit. Mach contours. 


Figure 13 2D Telescope inlet with one internal design. 
Numerical simulation results based on full Navier-Stokes 
ecuations using NASA CFL3d code. 

















